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We show how changes in unitarity-preserving boundary conditions allow continuous interpolation 
among the Hilbert spaces of quantum mechanics on topologically distinct manifolds. We present 
several examples, including a computation of entanglement entropy production. We discuss approx- 
imate realization of boundary conditions through appropriate interactions, thus suggesting a route 
to possible experimental realization. We give a theoretical application to quantization of singular 
Hamiltonians, and give tangible form to the "many worlds" interpretation of wave functions. 

PACS numbers: 



Quantum theory, starting with its name, has a repu- 
tation for introducing discreteness. It can, however, also 
blur it. Here we show how the behavior of quantum me- 
chanical models based on given manifolds can be interpo- 
lated continuously, as the manifolds undergo topological 
changes. 

The application of quantum mechanics to general rela- 
tivity suggests the possibility of violent fluctuations in ge- 
ometry under extreme circumstances, particularly in vir- 
tual processes at ultra-short distances. It has widely been 
speculated that changes in spacetime topology might oc- 
cur as a result. Unfortunately it appears very difficult 
to perform detailed or convincing dynamical calculations 
for quantum gravity, even in the possibly consistent form 
of string theory, or to adduce relevant phenomena. (The 
joining and splitting of strings themselves is a form of 
topology change, whose description involves the special, 
complicated apparatus of vertex functions and string field 
theory.) But as we will describe, in the relatively mun- 
dane context of quantum mechanics, topology change 
amounts to no more than a smooth change of boundary 
conditions in Hilbert space, or, equivalently, a change of 
parameters in the Hamiltonian. The concrete framework 
we present may suggest a less singular, more inherently 
quantum mechanical approach to topology change, and 
perhaps shed light on the emergence of spacetime itself. 
It is not implausible that systems of the sort we consider 
could be realized experimentally. 

The idea of changing topology by smoothly varying 
boundary conditions in Hilbert space was proposed in 
Ref. Ij and further developed in Ref. [2^. In [1 it was 
pointed out that topologically distinct boundary condi- 
tions consistent with unitarity correspond to different 
self-adjoint extensions [3 4 of the Hamiltonian, between 
which there exists (at least in simple cases) a continu- 
ous interpolation through the space of possible unitary 
boundary conditions. 

Here we take a concrete approach to topology change, 
in which boundary conditions are implemented in Hilbert 



space by projection operators, and develop simple exam- 
ples in detail. For the most basic example, the fission of 
one interval into two, we follow the behavior of the eigen- 
functions under adiabatic interpolation between the two 
topologies, and compute the entanglement entropy pro- 
duced by this process. We also consider an application 
to a singular Hamiltonian system, which admits multi- 
ple topologically inequivalent quantizations. Finally, we 
discuss the implementation of boundary conditions via 
direct modification of the Hamiltonian, and close with 
some comments regarding possible further applications. 

Two-Point Unitary Boundary Conditions: To begin, 
consider the quantum mechanics of a particle supported 
on a space parameterized as two rays — oo < x < a, 
—oo<x<b. We assume that the kinetic energy density 
is given, away from the endpoints, by jcJa^i/jp; a smooth 
potential is also allowed. The Schrodinger equation gives 
us a probability current j ~ i(il}* dx'^l^—dx^* ip)- To insure 
that no probability flows out at a, b, we need a linear 
boundary condition satisfying j(a) = j{b) = 0. Any two 
Robin boundary conditions 

= a4,{a) + P^ia) 

= ^m + S^{b) (1) 

with fixed real a, /3, 7, S do the job. Conventional Dirich- 
let (/3, (5 = 0) or Neumann (a, 7 = 0) boundary conditions 
are of course included as special cases. The conditions ^ 
also are sufhcient to give a good eigenvalue problem for 
the Hamiltonian, since they are appropriate to terminate 
a second-order differential equation. If we impose bound- 
ary conditions of the type ([!]) separately at both a and 
b, we will describe the dynamics of a particle that lives 
in two independent spaces, with a fixed (by initial condi- 
tions) probability to be found on either one. At another 
extreme, we can also obtain unitary dynamics and a good 
eigenvalue problem by imposing two "identification" con- 
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Indeed, from the point of view of the wave function 
Eqn. ([2]) ensures that a and b are in effect the same point 
(with reversed orientation for the second ray), so that we 
obtain the quantum mechanics of a hne. 

A more general framework includes both Eqn. ([T]) and 
Eqn. ([2]) as special cases, and allows smooth interpolation 
between then. To capture the underlying formal struc- 
ture, we define the vector 



V^(a) m §|(a) 11(6) 
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and matrix 
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With these notations, the current conservation condition 
reads 



uUu = 
vj^ Jv = 



(4) 



where in the second line we display the polarized form. 
Now we implement boundary conditions by projecting 
the allowed vectors v to be of the form ti = H^, with 
a fixed real projection operator 11 and an arbitrary nu- 
merical vector ^. In this formulation, for example, we 
represent Dirichlet boundary conditions and identifica- 
tion boundary conditions respectively with 
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We will have current conservation within the allowed 
Hilbert space if 

nVn = 







(7) 



To obtain a conventional eigenvalue problem we want 
two conditions, so that 11 will represent projection onto 
a two-dimensional subspace. Projection onto a subspace 
of larger dimension will not permit one to satisfy Eqn. ([7|. 
Projection onto a subspace of lower dimension leads to 
an unconventional, though not necessarily inconsistent, 
eigenvalue problem and Hilbert space. The antisym- 
metric matrix J implements a hermitean symplectic in- 
ner product, and maximal (here, two-dimensional) pro- 
jections satisfying Eqn. ([7| define Lagrangian subspaces 
with respect to J. 



Interpolation: We can interpolate continuously be- 
tween identification boundary conditions 11= 'At 6 — 
and Dirichlet boundary conditions II/j at = 
±7r/2, while satisfying the requirements of Eqn. ([t]) and 
rank(n) = 2 throughout, in the manner 
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(abbreviating cos 9 c, sin 61 — s). 11(6') implements 
projection onto the subspace spanned by (0,0,1,-1)^ 
and {c,c, s, s)'^ . Note that for intermediate values of 6, 
this subspace does not have a purely geometric interpre- 
tation in physical space. 

This basic joining - or, read in the opposite sense, split- 
ting - operation can be used as a module to implement 
more complex topological transformations, such as the 
pinching off of a closed "baby universe" - that is, a cir- 
cle - as smooth interpolations among consistent, reason- 
ably conventional quantum mechanical systems. Figure 
1 illustrates that geometric possibility. (A different in- 
terpolation between one circle and two appears in ^.) 





FIG. 1: Starting with a circle, breaking at P — )■ a, & and 
Q — ^ c,d, then joining a,c ^ R and b,d ^ S yields two 
circles. 



The concept of smooth interpolation between topolo- 
gies can be extended to higher dimensions. In that con- 
text, we need conditions on the wave function and its 
normal derivative. In Figure 2 we illustrate a continuous 
"surgery" leading from two spheres to a torus plus four 
waste disks. In these constructions, the interpolation an- 
gle 6 can, but need not, be a (spatially) constant function 
on each bounding circle. 

Groups and Geometry: We now briefly elaborate on 
the general framework wherein interpolation through 
boundary conditions is naturally discussed. (This struc- 
ture arises, and has been studied, in the context of quan- 
tum graphs [5]. Related structures appear in [1] and 0.) 

Consider a network of wires with junctions of arbitrary 
multiplicity. We want to insure that no net probability 
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FIG. 2: By passing from identification to Dirichlet bound- 
ary conditions on four small circles, discarding the interior 
disks, and passing from Dirichlet to identification conditions 
on two circular boundaries, we interpolate continuously be- 
tween quantum mechanics on two spheres and quantum me- 
chanics on a torus plus four disks. 



flows into a given junction where n segments terminate 
at points Xj. (Given p wires we have n = 2p possible 
terminals.) If we orient all derivatives to define inward 
normals (that is, oriented away from the termini), then 
the necessary and sufficient condition for current conser- 
vation takes the form 



w^Jv 
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where 



and similarly for w, with ipi — > 7/'2, and 
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where each entry denotes an n x n matrix. Again J 
implements a hermitean symplectic inner product. 

Now we want to insure Eqn. ([9| by making a projection 
n of the allowed vectors onto an n-dimensional subspace. 
So we require 



nVn = 



(12) 



If we have a solution 11, then we can look for others by 
varying the subspace, so 11 — )• UIIU~^. Without loss of 
generality, we can take 11 to be hermitean, and U unitary. 
Eqn. (12 1 goes over into 
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(13) 



If we write the [7 = 1 + cM for infinitesimal e then 



Eqn. ( 13 ) gives 



M = 



A Bn 
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where Bh and Ch are hermitean. Combining this with 
unitarity (i.e., M antihermitean) we find 



M = 



Aa Bh 
-Bh Aa 



(15) 



where Aa is anti-hermitean. This form of M spans two 
commuting U(n) Lie algebras, i.e. 



M+ = 



Aa ^iAa 
-fi-Aa Aa 



(16) 



Exponentiating, we see that the "unitary hermitean sym- 
plectic" group is isomorphic to U(n) x U(n), embedded 
in U(2n) as 
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We get the geometry of allowed subspaces as a coset 
space, by choosing a reference IIq and dividing out by the 
isotropy subgroup of U matrices for which U^Wi^U = Hq. 
Choosing pure Dirichlet 
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as reference, we find that the diagonal subgroup is the 
isotropy subgroup. Thus the geometry of the set of 
Lagrangian subspaces is the geometry of U(n) x U(n) 
divided by its diagonal subgroup, i.e. the coset space 
\]{n) X U(n)/U(n)A. For the projector we have 



V = [tpi{xi),'ipi{x2),...,ll^l{Xn),'>P'i{xi),^[{x2),...,lp'i{Xn)) 
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(19) 



where W = C/JC/a- 

If we restrict to boundary conditions that do not mix 
the wave functions with their derivatives, then we have 
_B = in Eqn. ( 15 ). This picks out the diagonal subgroup 
U(n)A of the general U(n) x U(n), with Ui = 112- This 
group is naturally associated with boundary conditions 
of the form 
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R 
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(20) 



with R, S projection operators that treat wave functions 
and derivatives separately. Returning to the general case, 
for current conservation we want RS = 0, and for a good 
eigenvalue problem we want rank(i?) + rank(>5') ~ n, so 
as projections R + S = 1. The diagonal subgroup will 
preserve the rank of R and S separately, so we have n -t- 1 
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distinct orbits. Geometrically, the orbits are Grassman- 
nians U(n)/U(fc) x U(n — k). 

For example, for n = 2 and rank(i?) = rank(S') = 1 we 
can take the reference 



Rq 
Sq 



1(11 

2 I 1 1 



W 1 -1 

2 I -1 1 



(21) 



which implements the identification boundary conditions 
of Eqn. For this choice of Ro,So the isotropy sub- 
group is generated by (e***, e*^) and {e^^'^^ , 6*^'^^), and the 
orbit geometry is SU(2)/U(1). For Rq the identity and 
^0 = 0, or vice versa, corresponding to pure Neumann or 
pure Dirchlet, respectively, we get an orbit consisting of 
a single point. 

Examples: For the purpose of studying dynami- 
cal topology change, we would like to consider time- 
dependent boundary conditions as providing a kind of 
external field. Below we will indicate how this can be 
done rigorously; here we provisionally assume adiabatic 
evolution, and check its internal consistency. Within the 
space defined by 11 the instantaneous Hamiltonian de- 
fines a normal eigenvalue problem, and the amplitudes on 
different levels obey the usual quantum dynamical equa- 
tions. 

Returning to the example of fission of the interval 
[0, 271"] into two equal intervals, with a free particle Hamil- 
tonian ^ = ~ ( ) J we can carry forward explicit calcu- 
lations in the adiabatic approximation. It is convenient 
to use coordinates < a; < 27r for the mother interval and 
0<j/ = x<7r, 0<2 = 27r — x<7ron the daughters. 
We maintain Dirichlet boundary conditions at a; = 0, 2tt, 
y = 0, z = throughout, and evolve from identifica- 
tion to Dirichlet atx = y = z — t:, following our earlier 
prescription. Eigenfunctions can be taken real, and then 
follow the forms (sin ky, =p sin kz). We have the boundary 
condition 

(^(y = tt), tp{z = TT),ip'{y = tt), i/j'iz = tt)) 
(X (a cos 6*, a COS 6*, a sin 6* -I- 6, a sin — fe) (22) 

For the upper sign (sin ky, — sin kz) the solutions are 
k ~ I where I is a positive integer, independent of 9, 
with a = 0. In terms of x these modes are simply sin kx. 
They are antisymmetric around x = tt, vanish there and 
are smooth. 

For the lower sign (sin ky, sin kz) the solutions have 
b = and fccot kn — tan 0. In terms of x these modes 
take the form 9{tt — x) sin kx + 9{x — tt) sin fc(27r — x). 
They are symmetric around x — t:, and have kinks there. 
At = we have k = 1 — 1/2, where I is a positive integer, 
and k evolves continuously until k = 1 aX 9 = —tt/2. 

Since the gap between states of the same symmetry 
never closes, adiabatic evolution of states is well defined, 
and will govern sufficiently slow evolution in the regu- 
lated dynamical framework described below. There is no 



geometric (Berry) phase, since we can use bases of real 
eigenfunctions throughout. 

The adiabatic evolution of the eigenfunctions has a 
simple intuitive interpretation in terms of node adjust- 
ment. Eigenfunctions of the first class (upper sign) al- 
ready have a node at the midpoint, so they are pre- 
adapted to Dirichlet boundary conditions, and do not 
change. Eigenfunctions of the second class (lower sign) 
have a local maximum or minimum at the midpoint. As 
they evolve adiabatically, that extremal value is pushed 
down toward zero. The interior nodes move slightly in 
the direction away from the midpoint, but do not disap- 
pear, and no additional nodes are created. 

Interesting complications occur if one considers fission 
away from the midpoint. Node induction implies that 
wave functions of both kinds will be adjusted. If the 
point of fission only deviates slightly from the midpoint, 
then among the low energy states, the modes with k = I 
are pushed into the narrower interval and their energy 
floats up slightly, whereas the modes with k — I — ^ are 
pushed into the wider interval, ending up with k that 
are smaller but close to I. The energy gap narrows as 
we approach the Dirichlet point, which requires slower 
evolution. In the other extreme, if we split at very small 
X = e, the low energy wave functions are easily adaptable: 
they simply move out of the interval [0, e]. 

Quantum Entropy and Entanglement: Our midpoint 
fission of the interval yields a mapping of Hilbert spaces 
L2([0,27r]) ^ L2([0,7r]) © L2([o,7r]). (Note: The direct 
sum appears here, not the tensor product.) In that no- 
tation, we have an informative basis 



1 / . . 2p - 1 , / 2 . 

(smpx + sm — - — x) — W — smpw®0 
2tt 2 ' V TT 



1 , 2p-l . 2 

Vp = . — (smpx — sm — - — x) — t- W — smpz 
V27r 2 V TT 

(23) 

with p running through the positive integers. (Note that 
since Up, Vp are not eigenfunctions until the fission is com- 
plete, their coefficients in a dynamical wave function may 
have complicated time dependence prior to then. The 
mapping of Eqn. ( 23 1 is defined structurally, rather than 



directly dynamically.) 

Now if we consider an observer confined to the first 
daughter interval L^([0,7r]), since all the possible opera- 
tions are also confined to it, the appropriate density ma- 
trix would have two separate sectors, corresponding to 
one-particle states sin py and the no-particle state. We 
can simply parameterize them with indices p, 0. Start- 
ing with the wave function J2p otpU^ 
density matrix 



3pVp we find the 



Piq-.P) 

P(o,o) 

^(9,0) 



p 

P{0,p) = 



(24) 
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This is not a pure-state density matrix, but it is close 
to one: a renormalized version of P(p,q), dropping the no- 
particle entries, would be a pure-state matrix. A nonzero 
quantum entropy is associated with it 

-Tr(plog2p) = - rlog2r - (1 - r)log2(l - r) (25) 

where r = X^pl/^pP- maximized at one bit when 

l/3pP = 5- This result makes intuitive sense, because 
von Neumann entropy is the intrinsic part of the uncer- 
tainty of a quantum state (as opposed the excess classi- 
cal entropy unveiled by measurement), and confinement 
to one subspace simply deprives the observer of the con- 
trol over the other intrinsically. More formally, implicitly 
underlying the knowledge of particle number is the fact 
that the observer has applied a measurement by Pi, P2, 
orthogonal projectors onto the two subintervals respec- 
tively, turning the original p to PipPi + P2pP2- Such 
measurements always increase entropy (if not leaving it 
unchanged) . 

Already in the two-particle case we find much richer 
structure. The mother Hilbert space is I/^([0, 27r]) (E) 
L2([0,27r]), which fissions into (^^([O, tt]) © ^^([o, tt])) ® 
(L2([0,7r]) © L2([0,7r])). There are four sectors in the 
product, corresponding to both particles in the first in- 
terval, the first particle in the first interval and the second 
in the second, and so forth. Using the Up, Vp basis we can 
write the initial state as 



r,s 

(26) 

When we form the density matrix for the first interval 
there will be four sectors, corresponding to both parti- 
cles in the interval, the first particle only, the second 
particle only, or no particle. Their coefficients will be 
quadratic in a,/3,7,5 respectively. The most interesting 
cases are when we have one particle, let's say the first. In 
an evident notation, we have for that sector the density 
submatrix 



(27) 



Since several values of I can make nontrivial contribu- 
tions, this will generally not be a (renormalized) pure- 
state density matrix, and we have true entanglement en- 
tropy. 

In this discussion we have assumed the two particles 
are distinguishable. Bosons or fermions can be handled 
without serious difficulty, by imposing symmetry or an- 
tisymmetry on a, l3, 7, S. 

Boundary Conditions As Limiting Interactions: 
Boundary conditions can be implemented formally 
as singular interactions. For example i^ia) = is 
enforced, for finite energy states, by the interaction 
AH — w|'0(a)p as w — > 00, corresponding to the 
potential V{x) = v6{x — a). The boundary condition 
Tpia) = tp{b) can be implemented through the energy 



term AH = v\ip{a) — ■ip{b)\'^, which corresponds to a 
nonlocal potential in the form 

(V|Ai/|V) = (28) 

V J dy J dx ^p* {y)(^S{x — a)5{y — a) + S{x — b)5{y ~ b) 

~S{x — a)6{y — b) — 6{x — b)S{y — a)^xp{x) 

or as an operator 

{AHi;){x) = 

v{6{x - a)V'(a) + 5{x - b)ij}{b) 

-5{x - a)i){b) - S{x - b)i,{a)) (29) 

Boundary conditions involving spatial derivatives bring 
in potentials with derivatives of delta functions. 
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FIG. 3: By allowing stubs, we can relax the boundary condi- 
tions at a, h, c, d while maintaining strict unitarity. The stubs 
need not be short, nor separate. 

For many reasons, including possible practical real- 
ization, it is interesting to consider relaxing the con- 
ditions on these singular interactions. Specifically, we 
would rather not insist on strict delta functions or infinite 
coefhcients. This is not unproblematic, however, since 
such compromises threaten unitarity. It is also desirable 
to have all the dynamics occur within a single master 
Hilbert space, rather than a time-dependent sequence of 
different function spaces. At some cost in elegance we can 
implement our constructions approximately, while main- 
taining unitarity exactly, in a single Hilbert space as fol- 
lows. Instead of using real terminal points, where we need 
strict boundary conditions, in our constructions, we let 
a, b, c, d, ... be interior points, each attached to a stub (See 
Figure 3). Instead of using 5 functions (or derivatives of 
5 functions) localized near these points in Eqns. (28 29) 



and their generalizations, we use regular approximations 
to those functions, e. g. Gaussians and their derivatives, 
with large though finite coefficients. Then wave functions 
that are supported in the region between the stubs will 
realize, approximately, the dynamics of topology change 
we proposed above. Since the boundary conditions are 
only enforced approximately some probability will leak 
into the stubs, but if the approximation to the singular 
interactions is accurate that leakage will be small. 

Application to Quantization: In a recent paper [6], 
two of us proposed a method for quantizing Lagrangians 
whose kinetic terms are of the form Lkin. = \x'^ — f 
Standard quantization procedures do not apply (for k, > 
0), since the associated Hamiltonian is singular, with 
cusps at p± = linking three smooth branches where 

—00 < p < p+, p- < p < p+, p^ < p < 00. (In the 
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interval between the cusps the Hamiltonian, as a func- 
tion of p, is three- valued.) After introducing the poten- 
tial V{x) = f and making the duality transformation 
X ^ p, p ^ —X, the system maps onto the description of 
a particle with conventional kinetic energy on a singular 
space. Introducing the variable 

X = 2; — + x^ for X ^ 

X = — a; + .T-|_ + .T_ for x- ^ x ^ x+ 

X = X + Xj^ — x^ for ^ X (30) 

allows us to unfold the space, so now the branches oc- 
cupy —00 < X < a;_, x_ < X < 2;+, a;+ < X < CO; B,ud 
the salient issue is what boundary conditions to impose. 
In [6], guided by unitarity and simplicity, we chose what 
we've here called identification boundary conditions, that 
is continuity of the wave function and its first derivative 
at x± . This leads to a consistent quantization of the orig- 
inal classical Lagrangian, wherein motion through the 
cusps, taking the particle from branch to branch, is unim- 
peded. 

On the other hand classical mechanics for systems of 
this kind suggests the interest of allowing "hard wall" re- 
versals of X at the cusps [7] , which in turn might suggest 
imposing il}(x- — e) = ip(x- -|-e) = — e) = '0(2;+ -I- e) 
at our four endpoints. These three conditions together 
with ^'{x--e)-ii'{x^+e)+^'{xj^-(i)-'\ij'{x+ + e) = 
yield a consistent, unitary theory ("ring quantization"). 
Another possibility, perhaps suggested by quantization 
using noncanonical variables |8], [6], is to impose Dirich- 
let boundary conditions il>{x- — e) = '0(x_ -I- e) = 
"0(2;+ — e) = "0(2^+ + f) = 0. Again, this is consistent 
with unitarity and realizes the classical correspondence 
principle away from the cusps. An alternative interpre- 
tation of the hard wall is that a particle starting on 
the branch between the cusps remains trapped there for- 
ever; this suggests the baby universe boundary conditions 
{il){x-+e) = ■0(2;+ ^e); ilj'[x^+e) = ■0'(2;+ — e)}, supple- 
mented with Dirichlet ip{x- — e) = "0(2;+ + e) = or iden- 
tification {■0(x_—e) = ■0(2;+-l-e); ip'{x-—e) — ip' {x++e)} 
conditions at the ends of the outer branches. 

These several choices - and others - all fit naturally 
within the general framework proposed above. They 
correspond to different choices of Lagrangian subspaces 
within the eight-dimensional space spanned by the values 
of the wave function and its derivative at the four end- 
points, with respect to the symplectic structure imposed 
by unitarity. Each can be reached from any of the oth- 
ers by interpolation through consistent quantum theories. 
The circumstance that a single classical theory can give 
rise to many possible quantum theories arises in QCD, 
where one has the 9 parameter [TD]; or, more simply, in 
the quantum mechanics of a particle on a circle, where 
one has the possibility of fractional angular momentum 

Comments: 

1. The central formal ideas of the preceding con- 
structions, including association of quantum the- 



ories with Lagrangian subspaces, smooth interpo- 
lation between theories based on topologically dis- 
tinct manifolds by motion in the consequent moduli 
spaces, and the realization of all these structures 
approximately through simple quadratic interac- 
tions, carry over to interacting many-body systems 
in a straightforward way. With some modification, 
they also apply to other differential operators, e.g. 
Dirac, Majorana, Yang-Mills, or Einstein (curva- 
ture). It would be interesting to elucidate the rela- 
tion between topology and the spectral fiow of such 
operators, including specifically index theory, from 
this perspective [12j . In a similar spirit, one could 
calculate charge flows [13]. 

2. Similar structures arise in quantum field theory by 
a slightly different route, as follows. In that context 
we consider wave functionals ^'(0(x)), where now 
the functions 0(2;) live on various manifolds with 
boundary. The Hamiltonian typically contains a 
kinetic piece which is formally 

ffgrad (X y"(V0)2 (31) 

operating on ^E*. To define an appropriate Hilbert 
space of functions for which i/giad is Hermitean 
by local (or nearly local) conditions on the bound- 
ary, we are led to impose boundary conditions of 
the kind discussed above. 

3. If we allow spatial variations in the surgery con- 
structions then we have, in view of our discus- 
sion in "groups and geometry" , the structure of 
a U(2) X U(2)/U(2)a gauged sigma model on the 
boundaries. Such sigma models and their gener- 
alizations will arise in compactifications of extra 
dimensions. With four noncompact dimensions, 
instanton-like configurations based on 7r3(U(n)) can 
mediate changes in the topology of the compactifi- 
cation manifold. 

4. The idea that quantum wave functions potentially 
describe many worlds gets realized in a precise and 
tangible form, as the quantum dynamics of (say) 
a single interval evolves into that of two separate 
intervals. Let us begin with a pure state on the 
mother interval. After the interval has fissioned, 
we still have a pure state on the two-interval mul- 
tiverse. But observers (or observables) confined to 
one interval will find it appropriate to use a mixed 
state, described by a density matrix that traces over 
dynamical variables localized on the other interval. 
An important point is that the dynamics of these 
mixed states will appear perfectly local and ordi- 
nary. In that sense the evolution has revealed that 
two distinct daughter worlds were implicit in the 
wave function of one mother. Note that the quan- 
tum entropy induced by the fission depends on the 
initial state. 
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